EPJ manuscript No. 

(will be inserted by the editor) 



Commensurate-incommensurate transition in two-coupled chains 
of nearly half-filled electrons 

M. Tsuchiizu 1 , P. Donohue 2 , Y. Suzumura ' , and T.Giamarchi • 

1 Department of Physics, Nagoya University, Nagoya 464-8602, Japan 

2 Laboratoire de Physique des Solides, CNRS-UMR 8502, Universite Paris-Sud, Bat. 510, 91405 Orsay, France 

3 CREST, Japan Science and Technology Corporation (JST), Japan 

4 Laboratoire de Physique Theorique, CNRS UMR 8549, Ecole Normale Superieure, 24 Rue Lhomond 75005 Paris, France 

Received: 

Abstract. We investigate the physical properties of two coupled chains of electrons, with a nearly half-filled 
band, as a function of the interchain hopping t± and the doping. We show that upon doping, the system 
undergoes a metal-insulator transition well described by a commensurate-incommensurate transition. By 
using bosonization and renormalization we determine the full phase diagram of the system, and the physical 
quantities such as the charge gap. In the commensurate phase two different regions, for which the interchain 
hopping is relevant and irrelevant exist, leading to a confinement-deconfinement crossover in this phase. 
A minimum of the charge gap is observed for values of tj_ close to this crossover. At large t± the region 
of the commensurate phase is enhanced, compared to a single chain. At the metal-insulator transition the 
Luttinger parameter takes the universal value K* = 1, in agreement with previous results on special limits 
of this model. 

PACS. 71.10.Hf Non-Fermi-liquid ground states, electron phase diagrams and phase transitions in model 
systems - 71.10.Pm Fermions in reduced dimensions (anyons, composite fermions, Luttinger liquid, etc.) 
- 71.30.+h Metal-insulator transitions and other electronic transitions 



1 Introduction 

Organic conductors Bechgaard salts, which consist of an 
array of one-dimensional chains, show various ordered states 
such as spin-Peierls state, spin density wave state and su- 
perconducting state at low temperatures [y . Besides these 
states of broken symmetry, the normal state above the 
transition temperature exhibits a remarkable electronic 
state associated with a charge gap. This charge gap is sup- 
posed to originate from the electronic interactions and the 
commensurate band-filling (1/4 filling) of the compound 
0. Indeed in one dimension, commensurate systems are 
Mott insulators. Such an interpretation of the charge gap 
has received support from recent optical experiments p, 
^| and transverse conductivity measurements 13]. In ad- 
dition, since these compounds are quasi-one-dimensional 
systems there is a competition between the one dimen- 
sional charge gap, that tends to localize the electrons, 
and the interchain hopping tending to make the system 
three (two) dimensional. This competition could lead to a 
confinement-deconfinement transition responsible for the 
difference of behavior between the TMTTF and the TMTSF 
salts [0,^] ■ Experimentally the confinement (deconfinement) 
is found in TMTTF (TMTSF) salts whose interchain hop- 
ping is smaller (larger) than the charge gap ||. From a 
theoretical point of view, the transition should take place 



when the interchain hopping renormalized by interactions 
is comparable with the single chain gap There is 

thus still considerable debate on how such a transition 
takes place, and what are the relevant energy scales. 

Unfortunately studying an infinite number of coupled 
chains is extremely difficult. It is thus worth to investi- 
gate these issues on a finite number of coupled chains, i.e. 
on ladder systems || . By studying explicitly two-coupled 
chains with a half-filled band, it has been shown that 
the interchain hopping becomes irrelevant, in the sense 
of renormalization group, for a charge gap larger than the 
interchain hopping It is reasonable to identify this 

relevance (irrelevance) of interchain hopping with the de- 
confinement (confinement) for an infinite number of chains 
since the interchain hopping is renormalized to zero in the 
limit of low energy for the irrelevant case. Of course there 
are some differences between the simplified ladder system 
and the infinite number of chains: (i) for the ladder the 
confinement-deconfinement transition is in fact a crossover 
as far as the ground state is concerned |TT1 , |T^| ; (ii) in the 
ladder both the relevant and irrelevant cases lead to an 
insulating state due to a gap in the total charge excita- 
tion. Thus in the ladder to obtain a metallic state, as 
observed in the experimental (infinite number of chains) 
compounds, it would be necessary to explicitly dope the 
system and go away from the commensurate filling. This 
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could be a way to mimic in this simplified model the small 
deviation of the commensurate filling due to warping of 
the Fermi surface perpendicular to chain direction ||. 

Even if it is known that upon doping the ladder be- 
comes metallic and the confined region is suppressed |l3| , 
it is yet very unclear what are the full properties of the sys- 
tem upon doping and what are the characteristics of such 
a metal-insulator transition. The purpose of the present 
paper is to examine these questions in the ladder sys- 
tem. We show that the metal-insulator transition in the 
ladder can be accurately described by a commensurate- 
incommensurate transition, and we determine its charac- 
teristics. The commensurate-incommensurate transition is 
well known in the classical case for one-dimension |T^ . p^| 
or quasi-one-dimensional [[l6| system. However the quan- 
tum case which corresponds to interacting electron sys- 
tems, is known only for one-dimensional case where 
its connections to the Mott transition have been investi- 
gated in details |l{|^0|,|]]. In the ladder, although both 
the half-filled commensurate system (spin ladder) and the 
extremely incommensurate case have been widely investi- 
gated gnnil 



comparatively little has been done 
on the metal-insulator transition close to half-filling. The 
ladder system with the umklapp scattering has been exam- 
ined in the chain basis, where the incommensurate phase 
shows no gap in the total charge fluctuation [Q. Stud- 
ies using either a mapping on a 50(8) symmetric model 
|^6|,^7|, onto a hard core boson system (28) or in the 
large interchain hopping limit |29|| show a drastic modi- 
fication of the universal properties of the metal-insulator 
transition compared to the single-chain case. The univer- 
sal properties close to half-filling have been checked nu- 
merically by DMRG |30|. We use here the bosonization 
technique and renormalization group to study the full 
problem as a function of the doping and the strength of 
the interchain hopping. This allows us to obtain the full 
phase diagram and in particular the interplay between the 
confinement-deconfinement at commensurate filling and 
the metal-insulator transition upon doping the ladder. 

The present paper is organized as follows. In Sec. II, 
formulation is given in terms of a bosonized phase Hamil- 
tonian and renormalization group equations are derived 
by assuming scaling invariance for response functions. In 
Sec. Ill, the phase diagram of the commensurate state and 
the incommensurate state is calculated by integrating the 
renormalization group equations. The charge gap is also 
estimated and the critical properties of the transition are 
given. In Sec. IV a summary and a discussion of the re- 
sults can be found. Technical details can be found in the 
Appendix. 



2 Formulation 

2.1 System at half-filling 

We consider two-coupled chains of a quarter-filled Hub- 
bard model with a dimerization given by 

« = -EE Et <+ c- 1 )^] ( c k* c i+^.< + h -' 

j (T=T4«=1>2 



-2*iEE ( c Li c ^>>2 + h.c.) 

3,1 



(1) 



where <r (=f, J, or +, — ) and I (= 1, 2) denote the spin and 
chain index, respectively, and td is the dimerization in the 
one-dimensional chains. After the diagonalization of the 
id-term, the kinetic term is written as 



k,(T,l 

-2^E[ d L,A^2 + h.c.] (2) 

k.a 

-2t± / ][yj. (T .i U k,<r,2 + b-C.] , 



where Sk = —2yt 2 cos 2 ka + t\ sin 2 ka |32) with the lattice 

constant a and dk.a-,1 ( u k,a,i) is the fermion operator for 
the lower (upper) band on the l-th chain. Note that the 
umklapp scattering is induced by the dimerization, which 
has an effect of reducing the quarter-filled band into an 
effectively half-filled one. Since we consider only the lower 
band, we use half-filling instead of quarter-filling in the 
present paper. By diagonalizing the interchain hopping 
term with the use of a^,a,C, = (— Cdfc,<r,i + dk,<r,2)/V2 (C = 
±), we obtain H Q = Efc^c e ( fc > 
energy dispersion is given by 



e(fc, C) = -2y i 2 cos 2 ka + t\ sin 2 ka - 2t ± ( . (3) 

Here we introduce a linear dispersion, e(k, C) — *■ vp(pk — 
&fc) where p is the index of the branch p = +(— ) for the 
right moving (left moving) electrons and the new Fermi 
point is given by = fcp — C2£_l/vf- We have neglected 
the independence of the Fermi velocity. After the bosoniza- 
tion of electrons around the new Fermi point fcpc, we in- 
troduce the phase variables defined by 



po i 



1 



where f(u,() = 1 (f°r v = p)> a (f° r v = c)j C (f° r 
v = C) and a( (for v = S), respectively. The phase vari- 
ables 9 P+ and 9 a+ (9c+ and 0s+), express fluctuations 
of the total (transverse) charge density and spin density. 
They satisfy the boson commutation relation given by 
[6 v +(x), 6 V '-(x')\ = iir8 v y sgn(x — x'). In terms of these 
phase variables, the field operator, ip p , a x 
ak,p,<r,C is given by pO 



1 



V2t 



exp (ipk FC x + iOp^x) exp (iir^, p ^ x ) . 



(5) 
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where a is of the order of the lattice constant and Pt aX = 
p/(2V2)[9 p+ + P 8 p _ + a(8 a+ +p6 a J) + C(0 C+ + P&C-) + 
<t((8s+ + p0s-)]- The phase factor, ir^, PtCr ^, in Eq. (§, 
is introduced to ensure the anticommutation relation for 
ip Pt<7 £ with different p, a and £. 

2.2 System close to half-filling 

In order to consider the system, which is slightly away 
from half-filling, we consider the following additional term, 



V2 



dx d x t 



p+ i 



(6) 



where \x is the chemical potential and fi = corresponds 
to the half-filling. We apply the transformation \/28 p+ — > 
V28 p + + qo% with qo = 4fj,K p /v Pl which leads to a mis- 
fit term q$x in the cosine term expressing the umklapp 
scattering. The phase Hamiltonian is written as, 



^ Air 

v=p,a,C,S 

9p 



-in 

dx 



dx 



-^{d8 v +f + K v {d6 v -f 



x 



cos • 



fJa 



An 2 a 2 



flu 



dx 



dx 



cos 





8* ± V 




x + 




VF J 


#3 + 


— cos 




8t± \ 


C+ 


x ~ 




v F ) 



cos 



V26» c _ 



V28r 



cos v / 2^s+ + cos V28g- 

V28 p+ + qax^j 
8ti 



cos( v / 26'c+ —a; 



cos 

cos V2&s + + cos v2#s- 



27r 2 a 2 



J dx cos V28 a+ 



8tj 



cos | v / 26'c+ a; 



cos V28c- 
cos V20S+ — cos \/20s- 



(7) 



where w p((T ) = Up[l +(-) U/irv?] 1 / 2 , v C (s) = «f, K p{a) = 
[!+{-) U/itvf}- 1 / 2 , K c(s) =l,g p = -g a = g ± = Ua and 
a coupling constant for the umklapp scattering is given by 
93 = -Ua(2t d /t)/[l + {t d /tf] §§. 

2.3 Renormalization group equations 

By utilizing a renormalization group method, we analyze 
Eq. (Q) where the nonlinear terms in Eq. (Q) are rewritten 



n 9 9 9vp-v'p> [ dx cos V2 6 vp cos \Pl 6v' p i , (8) 
2ir z a z J 

where 8 p+ = 8 p+ + q x/\/2, 8c+ = 8 p+ + 4\/2t±x/vF and 
9 up = 8 V p otherwise. The coupling constants are given 
by 9 P +,c± = T9 P +,s± = 9u, ±g<j+,c± = ~9a+s± = 9±, 
±9c±,s± = (g P + 5<t)/2 and ±gc±,S^ = ~{g P ~ 9*)/%, 
where each coupling constant is treated in the renormal- 
ization group method. The renormalization group equa- 
tions are derived from the response functions, which are 
given by (T T exp[i6 p+ (xi, n)] exp[~i8 p+ (x 2 , r 2 )]) By 
making use of the scaling of the cutoff (a — * a' = ae dl ), 
we obtain the equations as (Appendix A), 



dl 



t±=t ± - -G^^KcF^t^qoa) 



( G l+ : c+ + G c+,s+ + G c+,s-) K c Ji(8t±), 

(9) 



dl 



q a = q a - G 2 p+C+ K p F 1 (q Q a; 8t±) 

- (Gp+,c- + G>f,5+ + g2 p+s-) k p MQoa), 



dl K > = -~ K 



1 

2 



G 2 p+>c+ F (8t x ;q a) 



+ { G 1+,c- + G l+.s+ + G 'p+.s-) Jo(qou) 

1 



dl 



G 1+,c+ Ja{8t±) -\ 
~ G l+,s+ + G 1+,s- 



G l+.c- 



(10) 
(11) 
(12) 



■lY,[ K cM8t±)8p,+ -8p,-] (G 2 Cp 



s+ 



P =± 



~ v ' G Cp,S- + G a+,Cp) 

-\Y,[ K c F ^r,q a)Sp,+ 



p=± 



dl 



-Jo(q a)S p ^]G 2 p+ ^ Cp , 

K s = -\ Y ( K s s p,+ - S p-) g2 p+,s p Mqa a ) 



+ G l+,s P + G c+,s P Jo{&t±) 



G c-,s P 



(13) 



(14) 



dl 

d 
dl 



G u+,C± 
G u+,S± 



(2-K v - K%) G v+ , c ± 

—G v+ ,s+Gc±,s+ ~ G u+t s-Gc±,s-,(15) 

(2-K„-K±) G v+tS ± 

—Gv + ,c+ Gc+,s± Jo(8tj_) 

—G v +,c- Gc-,s± ! (16) 



d 



= f 2 - K* - K p Q ) G, 



*Cp,Sp' 



-G p+ ,cpG p+ ,sp> Jo(qaa) 
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-Ga+,CpGa+,Sp' 



(17) 



where v = p,a and p = ±. The quantities Fo(x;y) and 
Fi{x;y) are defined by 

Fo(x;y) = \[M\x + y\) + M\x-y\)} , 

Fi(x;y) = ^[Ji(\x + y\) sgn(x + y) 
+J 1 (\x-y\)sgn(x-y)] , 

and J„ is the n-th order Bessel function. In the above 
equations, the /-dependence is not written explicitly. The 
initial conditions are given by ^(0) = K v , G vp>v t p i(G) = 
9vp,v' P ' /2irvp, t±(0) = t±/W with W = vpa^ 1 , and <7o(0) = 
qo, respectively . The cutoff a can be related to the lattice 
constant by a — 2a/ir fjlj , which leads toW ~ nt/y/2 for 
small td/t. We take a=l. The renormalization group equa- 
tions for the velocity v v are discarded and the velocity v p 
and v a are set to vp. 

In the limit of t± — > these equations reduce to those 
of a single chain given by 

d jq a = q a-4 G 2 U J x {q a) , 



dl 
d G 



where K„ 



d 

Jl 
d 

Jl 
d_ 
dl 
1 



G u 

G a 

G x = 
— G n 



2G 2 u J (q a) 
2 G p G u , 



2G\ 



2 G a Gj 
and K„ 



(18) 
(19) 
(20) 
(21) 
(22) 



1 

-,s±, 



G a . In the above 
G± = ±G a+: c± — 
G„) = 



equations, G u = G p+ c± — TG 
-G a+iS ±, {G p + Gj) = ±2G c± ,s± and (G p 
=F2Gc±,s , =f- I* i s straightforward to derive Eqs. ([l8|)-(p2"| 
from the one-dimensional (ID) Hamiltonian given by (Ap 
pendix A) 



Hid = 



11. 

47T 



dx 



1 







1 




j dx 




47T 








dx d x 6+ 









K p (d x 9-Y 
+ K a {d x cf>. 



9u 



2w 2 a 2 

9± 
'2n 2 a 2 



dx cos20 + 



dx cos 20 + 



(23) 



where 8± and <f>± are the phase variables expressing the 
charge and spin fluctuations, respectively |Q . 

3 Commensurate and incommensurate states 

We examine the commensurate-incommensurate transi- 
tion by solving the renormalization group equations nu- 
merically. The dimerization is taken as td/t = 0.05 where 



o 



1 1 


1 1 1 1 TT T| 

(1)/ / 




/ 




// (3) \ " 




(4)\ \ 

i , i , \ i 



Fig. 1. The scaling flows of qo(l)a with fixed fi/t = 0.15(1), 
0.126(2), 0.125(3) and 0.1(4) for U/t = 5, t d /t = 0.05 and 
t±/t = 0.1. 



such a choice does not change qualitatively the results as 
seen later. The scaling quantity / is related to energy w 
and/or temperature by I = ln(W/u>) = \n(W/T). 

Equation (|lTJ) shows that the quantity qo(l) in the ab- 
sence of the interaction increases as qo(l) — qoe 1 , while the 
increase of qo(l) is suppressed by the presence of G p + : c± 
and G p +,s± coming from the umklapp scattering g u . In 
Fig. 1, qo(l) are shown with several choices of /i/i for 
U/t = 5, td/t = 0.05 and t±/t = 0.1. For fx/t = 0.15 
(curve (1)), the misfit term qo(l) increases rapidly and be- 
comes relevant. In this case, the umklapp scattering can be 
neglected since the Bessel functions Jo(qo<y) and J\(qoa) 
in Eqs. (S)-(17) become small due to large q^a. The quan- 
tity Kp(l) at the fixed point takes a finite value due to the 
relevant q$ot, and then the incommensurate state is ob- 
tained. For smaller value of the chemical potential given 
by fi/t = 0.1 (curve(4)), the quantity qo(l)a has a maxi- 
mum and reduces to zero with increasing I, showing the 
irrelevant qo(l)a. Such a behavior of qo(l)a leads to the 
commensurate state, which gives an insulating state due 
to the zero value of K p (l) in the limit of large We note 
that the commensurate state corresponds to that of a half- 
filled case. The commensurate-incommensurate transition 
occurs at the critical value given by fi/t = /x c /t(~ 0.125) 
(curve (3)). 



3.1 Half-filled case 

It is known that a transition from the relevant interchain 
hopping to the irrelevant one occurs with increasing t± 
at half- filling |],|l(]]. We show that such a result is also 
obtained for \i ^ 0, which still leads to the commensu- 
rate state. In Fig. 2, the quantity t±(l)/t is calculated for 
the commensurate state with fi/t = 0.1. With decreas- 
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/ 

Fig. 2. The scaling flows of t±(l)/t with fixed t±/t = 0.1(a), 
0.057(b), 0.056(c) and 0.02(d) for fx/t = 0.1, U/t = 5 and 
td/t = 0.05. 

ing the interchain hopping as t±/t = 0.1(a), 0.057(b), 
0.056(c) and 0.02(d), the hopping changes from a diver- 
gent behavior to a convergent one with zero value. One 
finds the relevant interchain hopping for large t± (curve 
(a)) but irrelevant one for small t± (curve (d)). The be- 
havior which separates these two cases is obtained at a 
critical value of t±_ = tx,c (curve (c)). For a relevant in- 
terchain hopping, one obtains an insulating phase with 
spin gap excitations which is called "D-Mott" phase |^G| 
with "C0S0" phase. The notation CnSm denotes a state 
with n massless charge mode and m massless spin mode. 
Such a spin gap is also known in the SO (8) Gross- Neveu 
model [ p6| , which includes an additional term given by 
9p+,a+ sinv^p-i- cosV^W- For an irrelevant interchain 
hopping, one can also expect a state with the spin gap 
due to the following fact. The phase has all charge ex- 
citations gapped, and because of such generated Heisen- 
berg exchange (J ~ t 2 ± /A p ) Jll| is equivalent to a spin 
ladder. Such a system is known to have all spin excita- 
tions gapped. We note that, for two-coupled chains, these 
states undergo a crossover since both states exhibit C0S0 
phase and there is no clear distinction between these two 
states |Il],[l2|]. However we may identify the irrelevant (rel- 
evant) t± as confinement (deconfinement) since such a re- 
sult could be expected for the case of many chains j}5| 
and higher dimension. 

3.2 Commensurate-incommensurate transition 

Based on Figs. 1 and 2, we examine quantities q and t±, 
which express effective quantities of qo and t± in the pres- 
ence of the interaction and the chemical potential. These 
quantities are estimated from qa — exp[— l q ] x 4 with 
q (l g )a = 4, and tf = texp[-h] with t±(h)/t = 1. One 



\il t 

Fig. 3. The effective quantity q vs. the chemical potential fi/t 
with fixed t±/t = 0, 0.05 and 0.2 for U/t = 5 and t d /t = 0.05. 
The thin dotted line denotes qo- The inset shows t Q f /t±, for 
t±/t — 0.05(dotted curve) and 0.20 (dashed curve). 



finds that q = qo and t c ® = t± for the nonintcracting case. 
In the present analysis, the effective quantity q is related to 
the carrier density n by n = q/ix (Appendix A). In Fig. 3, 
the //-dependence of q is shown with fixed t±/t = 0, 0.05 
and 0.2 for U/t = 5 and td/t = 0.05. The quantity q, which 
is suppressed by the umklapp scattering, becomes zero 
in the commensurate phase. The thin dotted line denotes 
q(= qo) in the absence of the umklapp scattering. The crit- 
ical value fi c for the commensurate-incommensurate tran- 
sition is not monotonical as a function of t± since \x c for 
t±/t = is larger (smaller) than that for tj_/t = 0.05 (0.2). 
The explicit ^-dependence of \i c is evaluated in Fig. 4. In 
the inset of Fig. 3, the effective interchain hopping nor- 
malized by t± is shown for t±/t = 0.05 (dotted curve) and 
0.2 (solid curve) where t Q f /t± is symmetric with respect 
to fi = 0. For t±/t = 0.05, t ^ /t± becomes zero for small 
n/t( ^ 0.11) and exhibits the irrelevant interchain hopping 
while the interchain hopping for t±/t = 0.2 shows always 
relevant one. 

From the results of Fig. 3, we examine the bound- 
ary between the commensurate state (Ci and Cn) and 
the incommensurate state (IC) and also the boundary be- 
tween the relevant interchain hopping and irrelevant one. 
In Fig. 4, the phase diagram of these states is shown for 
U/t = 5 and t^/t = 0.05. The solid curve denotes the 
boundary between the incommensurate state and the com- 
mensurate state. The incommensurate state corresponds 
to the metallic state. The commensurate state is the (Mott- 
)insulating state where the dashed curve in the commen- 
surate state denotes the boundary between the relevant 
interchain hopping ( "D-Mott phase" ) |ll| and the irrele- 
vant one ( "Confined phase" ) . The critical value shown by 
the arrow on the t^-axis is given by tj_ Jt ~ 0.072 [|lj 



6 



M. Tsuchiizu et al.: Commensurate-incommensurate transition in two chains of nearly 1/2-filled 




Fig. 4. The phase diagram of commensurate (Ci and Cu) 
and incommensurate (IC) states on the plane of fi/t and tj_/t 
for U/t — 5 and td/t = 0.05. The region Ci (Cu) denotes 
"Confined phase" ("D-Mott phase"), i.e., the commensurate 
state with the irrelevant (relevant) interchain hopping. The 
region IC corresponds to "Luther-Emery liquid". The arrow 
on the horizontal (vertical) axis shows the critical value fie/t 
(tx c /t) separating the respective regions. 




Fig. 5. The fj_-dependence of the charge gap A p /t for fixed 
n/t = 0, 0.1 and 0.14 for U/t = 5 and t d /t = 0.05. 



acquire a mean value (with opposite signs) so that the 
umklapp term reduces to: 



9u sin 



V28 



cos 



V29r 



— (cos 



y/20. 



s+, 



(24) 



|35| . The critical value shown by the arrow on the /^-axis 
is given by n®/t ~ 0.154, at which the solid curve merges 
with the dashed curve. The interchain hopping is always 
relevant in the incommensurate phase since the umklapp 
scattering becomes irrelevant in the limit of low energy 
as seen from Eq. (Q). The incommensurate phase shows 
no gap in the total charge fluctuation due to the relevant 
/i p2| , while the gaps still exist for other fluctuations due 
to the relevant t± . The incommensurate state corresponds 
to a "C1S0" d-wave superconductivity |^5| and is called 
"Luther-Emery liquid" [^9| . The solid curve corresponding 
to the critical value fj, c shows non-monotonical behavior as 
a function of t±. With increasing t±, fi c has a minimum 
at t±/t ~ 0.07 and becomes larger than yP c . The decrease 
of fi c for small t±/t( ^ 0.07) originates in the fact that the 
effect of the umklapp scattering becomes weakened for fi- 
nite t± due to the misfit (8t±/v-p)x in <? u -term of Eq. (Q). 

For large t±_ and small g u (<?„ is controlled by id which 
is small here) we can explain the behavior by a qual- 
itative analysis of the renormalization group equations. 
The strong interchain hopping between the two interact- 
ing chains opens a gap in all sectors except for the total 
charge sector. We may define a scale l\ where the cou- 
plings in the a, C, S sectors have reached a value of order 
one. This scale will depend on t± and the value of the 
bare interaction but not on g u . Up to h, g u will renor- 
malize by some finite multiplicative constant that will not 
affect the asymptotic dependence of the charge gap. How- 
ever above l\ the scaling dimension of g u is 2 — K p , instead 
of 2 - 2K p if t± = 0. Indeed cos{V26 s +) and cos(V26» c _) 



As a first consequence this means that the power law de- 
pendence of the gap for asymptotically small umklapp 
scattering is enhanced by interchain hopping from A 1 ^ oc 
^1/(2 2K P ) tQ ^ ^ gi/(2 k p ) ^ second physical con- 
sequence is for the commensurate-incommensurate tran- 
sition. In the absence of chemical potential of Eq. (|J) , the 
operator g u in Eq. (|24|) would be relevant when K^ =0 < 
2. As is standard for the commensurate-incommensurate 
transition ^,|), the addition of a chemical potential of 
Eq. (||) leads to a universal value of the Luttinger ex- 
ponent close to zero doping of K* = K£ =0 /2 — 1. It is 
also easy to see that the charge excitations connecting two 
minima of the potential, Eq. ( p4[ ) , correspond to a charge 
+2. We thus recover quite generally the results that were 
established in the various previous limits p7|, ^ , |29| ,p0|. 

Next, we estimate the magnitude of the charge gap 
on the plane of /i/t and tj_/t of Fig. 4. In a way similar 
to the previous calculations in the case for /i — P,[io|], 
the charge gap A p is calculated from the renormalization 
group flow where A p /t = W exp[— I a] and I a is defined 
by K p (l A ) = 1/2. We set A p = for K p (oo) > 1/2. In 
Fig. 5, the t^-dependence of A p is shown for fi/t = 0, 0.1 
and 0.14. When t± increases, the charge gap with /i = 
decreases and has a minimum at t±/t = 0.154. A similar 
dependence is found in the independence of fi c in Fig. 4 
where the location for the minimum is the same within the 
numerical accuracy. The identification of the charge gap 
for the commensurate case A p (fi = 0) with /i c can be ex- 
pected for the commensurate-incommensurate transition 
on general grounds. Indeed the charge gap is the smallest 
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U/t 

Fig. 6. The [//^-dependence of (solid curve) and t*± c (the 
dashed curve) for td/t = 0.05, which are shown by the arrows of 
the horizontal axis and the vertical axis in Fig. 4, respectively. 
The inset shows the corresponding id-dependence for U/t = 5. 



energy one must pay to inject an excitation which car- 
ries a charge. When the chemical potential reaches the 
charge gap, charged excitations are injected in the sys- 
tem, the system is incommensurate. The non-monotonical 
dependences of A p on t± indicate a crossover from the 
irrelevant interchain hopping to the relevant one with in- 
creasing t±. For [ijt = 0.1, A p is suppressed but is still 
similar to that of fi = 0. However the independence of 
A p for fi/t = 0.14 is qualitatively different from others. 
The metallic state with A p /t = appears in the interval 
region of 0.028 < t±/t < 0.2. This region corresponds to 
the incommensurate state in Fig. 4. We note that numer- 
ical integration of the renormalization group flow shows 
that Ap as a function of /i decreases monotonically to 
zero at the commensurate-incommensurate transition, in 
agreement with the annihilation of the gap expected for a 
commensurate- incommensurate transition. 

Finally we examine the U and id dependences of the 
critical values, /i° and v\_ c , which are shown by the arrows 
in Fig. 4. In Fig. 6 the [/-dependences of /i" and t_\_, c ar e 
shown by the solid curve and the dashed curve respectively 
for td/t = 0.05. In the inset, the ^-dependences of the 
corresponding quantities are shown by solid and dashed 
curves for U/t = 5. These U and id-dependences show that 
/Lt°Aj_ c — 2.2 for the present choice of parameters. Thus a 
rescaled phase diagram, which is independent of U and t^, 
can be obtained from Fig. 4 by using the rescaled variables 
and tx/*x c- Our renormalization group procedure 
correctly gives a value for /i" identical, within the scale 
of Fig. 4, with A p u . The quantity A p D denotes A p of the 

single chain where A p D ~ W(g u /W) 1 ^ 2 ^ 2KpS> for small 
9u @- We note that the boundary between a relevant t± 
and an irrelevant t± is determined more accurately by the 



competition between t e ^' and Ap^ where t°^'° denotes 
the effective interchain hopping energy for = ji = 
36[|35|l . The quantity t e F'° is given analytically by t e ±'° = 



t A _(t A _/t) a °/^ a °'> with a 



-4)/4 



4 Summary and discussion 

In the present paper, we investigated the properties of 
two-coupled chains with interchain electron hopping and 
a filling close to half-filling. By using bosonization and a 
renormalization group method we obtained the full phase 
diagram of the system. There is a metal-insulator transi- 
tion for a critical value fi c of the chemical potential, de- 
scribable by a commensurate-incommensurate phase tran- 
sition on the bosonized Hamiltonian. The critical value 
fi c , shown in Fig. 4, exhibits a non-monotonical depen- 
dence on the perpendicular hopping ij_. The minimum of 
\i c occurs for values of t± close to the ones for which a 
confinement deconfinement crossover takes place for the 
commensurate case. For large t± the relevance of inter- 
chain hopping reinforces the commensurate character of 
the system leading to an enhanced commensurability gap. 
The crossover line separating the regions of irrelevant and 
relevant interchain hopping (confinement-deconfincmcnt 
line) merges with the boundary between the commensu- 
rate state and the incommensurate state at // = /i", the 
critical chemical potential for a single chain (t± = 0). We 
found that the phase diagram of Fig. 4 becomes almost in- 
dependent of parameters such as interactions and dimer- 
ization when expressed in terms of the rescaled variables 
n/Hc (t±/t° ± ). In addition, given the form for the Hamil- 
tonian, we could show that at the limit of small doping, 
the Luttinger liquid parameter takes the universal value 
K* = 1, thereby confirming the results g^'U'll'Sl 
obtained on specific limits of the model. 

Finally let us discuss the interchain exchange inter- 
actions in the commensurate confined phase (irrelevant 
interchain hopping). Even when irrelevant, the interchain 
hopping generates two particles and particle hole hopping 
0. Within the present formalism, the two particle inter- 
chain hopping can be taken into account by starting from 
the chain. The Hamiltonian corresponding to the umklapp 
scattering, H u , in Eq. (0) is rewritten as 



P .i J 

P J 

p J 

+9u^2 / dx ^,T4^p,l,i ^-p>-1> 2 ^-p.T,2 + h.c. 



(25) 
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where ip p ,a,i = L~ x / 2 ^ fc e lkx dk, P ,a,i with I = 1,2 denoting 
the chain index and p = +(— ) corresponding to the right 
(left) moving electrons. The coupling constants are given 

by gi D = (g P +,c+ + g P +,c- - g P +,s+ + g p +,s-)/4, J uz = 
{g P +,c+-g P +,c--g P +,s+-g P +,s-)/^, J u ± = -(g p +,c+- 
g P +,c-+g P +,s++g P +:S-)/4&ndgP t = (g p+tC ++g p +,c-+ 

g P +,s+ - 5 P +,s-)/4, where gl D and J uz (J u ±) denote the 
intrachain umklapp scattering and the interchain umk- 
lapp exchange, and denotes the pair tunneling between 
chains. Other terms in the Hamiltonian can be rewrit- 
ten in a similar form. The initial values of the interchain 
exchange and that of pair tunneling are zero since the 
initial values of the umklapp scattering in the renormal- 
ization group equations are given by g p +.c+ = g P +.c- — 
~g P +,s+ = g P +.s- = gu- Both of these interactions are 
generated by the renormalization. 

It would be a interesting problem to investigate, by 
taking the higher order terms in the renormalization group, 
the consequences of these terms on the full phase diagram 
investigated in this paper. 
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the Ministry of Education, Science, Sports and Culture (Grant 
No.09640429), Japan. 



where Z = Tr exp(— J dr 7Yid) and G u = g u /(2itv p ). In 
Eq. (|28|), the new phase variable 9+(= 9+ — 2ttK p Xx/v p ) 

has been introduced and rewritten as 9+ — ► 9+. Then 
Eq. @ leads 

1 4 T f 

n = —q H ~G u K p — / dxdrx (sin (29+ + q x)) . 

2tt is or L J 

(29) 



Here it is worthwhile noting that Eq. (29) is compared 
with 



M = 



2K n 



-^ttG u — [ dx dr x (sin (26+ + 2-Knx)) . 
a 2 L J 



(30) 



which is obtained by using the Legendre transformation, 
i.e., by calculating the value of the chemical potential at 
fixed carrier density n p3|,^9|. 

After a straightforward calculation of Eq. (p9|), one 
finds, up to the lowest order of G Ul 



n = 7Tl° G u K p 



dr / r \ 2 ~ iK P , . , , 
- - J x (q r). 31 

a Va/ 



By assuming that the quantity n is simply scaled as n(l) = 
n e l [ ^3| with the transform a — > a e l , the renormalization 
group equation for q is obtained as, 



A Renormalization group equations 

In this section, we derive the renormalization group equa- 
tion for /i. 

First, we treat the system with a single chain where 
the phase Hamiltonian is given by 



n 



ID 



47T 



dx 



1 



(d x e+f 



gu 



2n 2 a 2 



dx cos (26*_| 



K p (d x 6-) 
\-q x) . 



(26) 



The quantities v p and K p are the same as two-coupled 
chains, and [9+(x), 6-(x')\ = Z7rsgn(x — x') Q and qo = 
(AK p /v p )n. Here the expectation value of the carrier den- 
sity, n, can be evaluated as 



2K n 



-M 



1 T 



dx dr (d x i 



(27) 



where the factor 2K p /irv p in the first term of r.h.s. corre- 
sponds to the compressibility in the absence of the umk- 
lapp scattering. The second term of Eq. (^7|) can be cal- 
culated as follows, 



dx dr (d x 9^ 



1 



:Tr 



Z 
4G U 



d_ 



exp 



dr Tim + A / dxdr d x 



K p J dxdr x (sin (26^ 



qox)) 



A=oJ 
(28) 



— q a = q a - AG\ K p Ji(qaa) 



(32) 



The renormalization group equations for K p and G u can 
be obtained in a way similar to Ref. 119] as 



-K p = -2GlK 2 p J (q a) 
~JjGu = (2 — 2K p ) G u . 



(33) 
(34) 



By integrating this renormalization group equation, the ef- 
fective quantity of go can be estimated from qa = ccxp(— l q ) 
with q (l q )a — c, where c is numerical constant of the or- 
der of unity. One finds that the quantity q can be related 
to carrier density n by n — q/(2w) from Eq. (^l|). 

Next we calculate the case of two-coupled chains. The 
expectation value of the carrier density, n, can be evalu- 
ated as 



AK P V2T 



v 2 T f 

A* H t dxdr (d x 9 p 

ir L J 



(35) 



where the first term of r.h.s. becomes twice as that in 
Eq. (|27|), since here we consider two chains. From a pro- 
cedure similar to the single chain, Eq. (|29|) is replaced by 



" = T'/n - ■^G p+tC + K p ^ / d.rdr .r 



1 -I 
— ( 

7T 7ra^ 



{ sin I V29 p + + qax) cos ( V29c+ x 



8fj 
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4 T 
— nGp+.c- K p — I dxdrx 



x y V29 p+ + q a xj cos \f2Q C - 

4 T f 

H ^G p+t s+ K P T / dxdrx 

■na z L J 



4 
■na 



(V29 
T 



i sin yyiu p+ + q x j <ios V29s+ 



jG p+ , S - K p — I dxdrx 



sin ( V%0 p+ + qox ) cos V29 



(36) 



The coupling constants that appear in Eq. ( p6| ) are those 
including the misfit qox in the cosine potential of Eq. ([?]). 
The calculation in the lowest order of perturbation yields, 



7T na 



dr 

a \a 



(D 



Fi (g r; 8i_i_r/u F ) 



f 



7ra 
f 



^ / r \ 2-2K„-2/K c 



dr / r \ 2—2K p —2Ks 



-/!(£) 



^ r /r v 2-2K p -2/K s 



Jx(qor) 
Ji(qor) 



(37) 



where Fi(x;y) = [Ji(\x+y\)sgn(x+y) + Ji(\x-y\)sga(x- 
y)]/2. The infinitesimal transform of the cutoff a — > a' — 
a e dl in Eq. ( jfi7\) leads the renormalization equationjEq. ( |l0| ) 
The renormalization group equation for Eq. (H), can 
be obtained in a similar way. The equations for the other 
coupling constants, Eqs. (|lll)-(|l7|) are also obtained in a 
way similar to Ref. p0|. 
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